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Abstract 

We theoretically analyze the influence of chemical disorder on magnetic anisotropy in Gei_ x Mn x Te semicon¬ 
ductor layers known to exhibit carrier-induced ferromagnetism and ferroelectric distortion of rhombohedral 
crystal lattice. Using DFT method we determine the local changes in the crystal structure due to Mn ions 
substitution for Ge and due to the presence in Gei_ x Mn x Te of very high concentration of cation vacancies. 
We calculate the effect of this structural and chemical disorder on single ion magnetic anisotropy mechanism 
and show that its contribution is order of magnitude smaller as compared to magnetic anisotropy mechanism 
originating from the spin polarization induced by Mn ions into neighboring Te and Ge ions. We also discuss 
magnetic anisotropy effects due to pairs of Mn ions differently allocated in the lattice. The spatial averaging 
over chemical disorder strongly reduces the strength of this magnetic anisotropy mechanism and restores 
the global rhombohedral symmetry of magnetic system. 
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1. Introduction 

GeTe belongs to the IV-VI family of narrow gap 
semiconductors. At temperatures below To = 670 K 
its crystal structure has rhombohedral symmetry 
Cs v which can be characterized by three parame¬ 
ters: the lattice constant ao ~ 6A, the rhombo¬ 
hedral angle a ~ 88.3° and the relative displace¬ 
ment of Te and Ge sublattices along [111] crystallo¬ 
graphic direction. This displacement which is equal 
to aoV3 t, where r ss 0.03, causes GeTe to be fer¬ 
roelectric. 

Manganese ions Mn 2+ introduce in Gei_ x Mn x Te 
local magnetic moments and this IV-VI diluted 
magnetic (semimagnetic) semiconductor (DMS) ex¬ 
hibits ferromagnetic transition with the Curie tem¬ 
perature Tc(x) up to 190 K 0 1- Due to si¬ 
multaneous ferroelectric and ferromagnetic order¬ 
ing Gei_ x Mn x Te belongs to the class of multifer- 
roic materials, which may become very interesting 
for possible future applications. 

Our previous theoretical analysis of the phys¬ 
ical mechanisms of magnetic anisotropy in 
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Gei_ x Mn x Te [3|, hereafter referred as paper I, 
was inspired by recent, unexpected results of 
magnetization and ferromagnetic resonance exper¬ 
iments which clearly showed that the easy axis 
of magnetization in thin layers of Gei_ x Mn x Te, 
for the manganese content of about 10 at. %, is 
perpendicular to the layer S 0, I, i 0, i. In 
ferromagnetic thin layers the shape anisotropy 
usually dominates with easy magnetization axis 
located in the layer plane. This is indeed ob¬ 
served in Gei_ x Mn x Te layers for higher content 
of manganese, above 20 at. %. For these layers 
the X-ray diffraction analysis (XRD) shows that 
the increase of Mn content results in transition 
from rhombohedral to cubic (rock-salt) crystal 
structure. 

In the paper I we performed extensive ab 
initio density functional theory (DFT) calcula¬ 
tions of magnetic anisotropy energy (MAE) in 
Gei_ x Mn x Te. The aim was to elucidate micro¬ 
scopic mechanisms responsible for the magnetic 
anisotropy. We concluded that the main contri¬ 
bution to MAE is given by the spin polarization 
induced in Te and Ge neighbors of Mn, and not 
by the Mn ions themselves. The Mn spin polarizes 
its neighborhood, the spins of Ge ions are approx- 
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imately parallel and the spins of Te ions approxi¬ 
mately antiparallel to the direction of the Mn spin. 
The degree and the spatial range of spin polariza¬ 
tion is relatively small in insulating crystal but it 
strongly grows with the increasing concentration of 
holes. Because of spin-orbit (S-O) interaction, the 
directions of Te and Ge spins forced by the direc¬ 
tion of Mn spin do not, in general, correspond to 
minimum of the total energy of the system. This 
forcing of Ge and Te polarizations by the Mn spin is 
the principal reason why variations of the directions 
of Mn spins result in changes of the total energy. 
Apart from the identification of the main physical 
mechanism responsible for magnetic anisotropy we 
also showed that MAE is determined by the hole 
concentration, macroscopic crystal structure, and 
the local configuration of Mn. 

The paper I was devoted almost exclusively to the 
situations where there was only one Mn ion in 2 x 
2x2 rhombohedral supercell containing 64 atoms. 
In all cases, the angular dependence of MAE on the 
Mn spin direction characterized by the angle 9 was 
well approximated by the formula 

Ea{9) = a, 2 c (cos20 — 1). (1) 

Here 9 is the angle measured in the (-110) plane 
from the [111] direction. Thus, for 0 = 0 and 9 = 
90°, the Mn spin is parallel to [111] and to the [11-2] 
crystallographic directions, respectively. The same 
definition of 9 is used in the present paper. 

One should keep in mind that a supercell contain¬ 
ing one Mn ion simulates perfect, translationally 
invariant crystal with rhombohedral symmetry. In 
actual random substitutional alloys, this symmetry 
is locally destroyed. 

Although macroscopic structural, mechanical, 
electrical or magnetic properties can have the global 
rhombohedral symmetry as confirmed by XRD, at 
the unit cell scale this symmetry is destroyed. This 
is due to the randomness of the occupancies of lat¬ 
tice sites in the chemically mixed cation sublattice 
(which is sometimes called chemical disorder), and 
to the presence of native defects, in our case the 
cation vacancies. The random positions of man¬ 
ganese ions together with different ionic radii of Ge 
and Mn lead to local deformation of the lattice. The 
local deformations are also caused by the germa¬ 
nium vacancies. Each cation vacancy, like in other 
IV-VI semiconductors, delivers two holes, thus in 
real sample with the hole concentration of the or¬ 
der 10 21 cm' 3 we expect 5 x 10 20 cm -3 of vacancies. 


For comparison, the Mn concentration x = 0.01 cor¬ 
responds to about 1.9 x 10 2 ° of Mn ions per cubic 
centimeter. 

The influence of the chemical disorder on the 
magnetic anisotropy in Gei_ x Mn x Te is the topic 
of the present paper. The paper is organized as 
follows. 

In Section [2] we analyze the single ion magnetic 
anisotropy (SIMA) in Gei-^Mn^Te. In the paper I 
we explicitly showed that SIMA in the conducting 
Gei_ x Mn x Te, is of minor importance compared to 
the anisotropy caused by spin polarized free car¬ 
riers. However, the calculations were performed 
for systems with manganese ion’s neighbourhood 
of perfect rhombohedral symmetry. On the other 
hand it is well known from numerous electron para¬ 
magnetic resonance experiments (EPR) that the 
ground state splittings and consequently magnetic 
anisotropy properties of Mn 2+ ion strongly de¬ 
pend on the symmetry of its surrounding. Previ¬ 
ous calculations performed for Pbi_ x Mn x Te and 
Sni-z-yPhyMn^Te show that the ground state 
splitting of Mn 2+ ion in disordered environment 
may be of the order of 1 Kelvin. Because the dif¬ 
ference in ionic radii of Ge and Mn is much smaller 
than that of Pb and Mn, the local lattice distortions 
(which strongly influence the magnetic properties of 
Mn) are also smaller. Thus, we do not expect that 
SIMA driven by the microscopic lattice distortions 
are dominant in Gei_ x Mn x Te. On the other hand, 
the band structure of GeTe is significantly differ¬ 
ent from those of other IV-VI semiconductors due 
to ferroelectric properties of the crystal. That is 
why direct calculations of SIMA in Gei_ x Mn x Te 
are necessary. 

In Section [2] we briefly describe mechanism re¬ 
sponsible for SIMA and the method of calculations. 
The obtained results confirm that even in the pres¬ 
ence of disorder, in conducting Gei_ x Mn x Te SIMA 
is of secondary importance. 

Section [3] presents results of ab initio calculations 
of MAE of pairs of Mn ions in Gei_ x Mn x Te. All 
nonequivalent positions of two Mn in the 2x2x2 
supercell are taken into account. We also show a 
few results for three Mn ions in the supercell. In 
the calculations the spins of Mn ions are always 
parallel. This is justified, because according to the 
results of Section II the energies connected with lo¬ 
cal magnetic anisotropy easy axes are small, and in 
the first approximation the possible noncollinear- 
ity of Mn spins can be neglected. We explicitly 
show that the direction of easy axis of magnetiza- 


2 


tion strongly depends on the relative position of Mn 
ions, and MAE is not described by Eq. ([1]). For cho¬ 
sen two configurations of two manganese ions we 
analyze contributions of kinetic, electrostatic and 
exchange-correlation energies to MAE. It is shown 
that the Mn ions introduced to GeTe causes signifi¬ 
cant perturbation in spatial distribution of electron 
charge and, consequently, significant change of elec¬ 
trostatic potential which, via the S-0 interaction, 
influences MAE. 

In Section U we summarize the main conclusions 
of the paper. 

2. Single ion magnetic anisotropy 

2.1. Theory 

According to the Hund’s rule, the ground state of 
Mn 2+ ion is 6 S, with the orbital momentum L = 0 
and spin S = 5/2. As the orbital singlet, it does 
not interact with the crystal environment, and its 
ground state without external magnetic field is six¬ 
fold degenerate. This model of Mn 2+ ion explains a 
number of phenomena, e.g., the results of magneti¬ 
zation measurements in diluted magnetic semicon¬ 
ductors. However, if we consider problems in which 
very small excitation energies are of importance, 
this model is no longer valid. An example of such a 
problem is the splitting of the ground state of Mn 
ion in a crystal held measured in EPR experiments. 
In the theoretical analysis for zero magnetic held, 
the Mn 2+ ion is described by a 6 x 6 matrix, Hmm ', 
which plays a role of effective spin Hamiltonian Q. 
The subscripts —5/2 < M,M' < 5/2 denote the 
projections of Mn spin on the quantization axis. 
For example, for a perfect cubic symmetry of Mn 
surrounding 

H=^St + S* + Si). (2) 

The lowest excitation energies of the ion (i.e., the 
eigenvalues of Hmm >) strongly depend on the sym¬ 
metry of its neighborhood, and are in the range 
10 -4 — 10 _1 K. If the Hamiltonian Hmm 1 is not di¬ 
agonal, it leads not only to the ground state split¬ 
ting but also to the magnetic anisotropy of the ion 

itio]. 

Density functional theory, which is one of the basic 
theoretical methods used to study condensed mat¬ 
ter properties, is not suitable for this problem be¬ 
cause DFT calculations determine the electron den¬ 
sity in the ground state, and provide no direct in¬ 
formation about the excited states of the system, 


particularly in situations where the energies of the 
excited states are very small. In this case one should 
use methods based on perturbation theory. 

There are two mechanisms known in the litera¬ 
ture which lead to the ground state splitting: in¬ 
teraction with excited states of 3 d 5 configuration, 
and hybridization between 3 d and orbitals of the 
surrounding ions. In the first mechanism Pj) the 
internal manganese S-0 interaction couples e S with 
the excited 4 P states, due to which the ground state 
of Mn is no longer a pure S state with total angular 
momentum L = 0. Consequently, the interaction 
with the crystal environment becomes possible. As 
a result, ground state splitting and nondiagonal ef¬ 
fective Hamiltonian Hmm 1 are finite. 

In DMS, the second mechanism, i.e., the hy¬ 
bridization between 3d(Mn) orbitals and the or¬ 
bitals of the neighboring ions, plays a more im¬ 
portant role. In this case, a system consisting of 
3d(Mn) electrons and valence band carriers is con¬ 
sidered. In the ground state of this system, which is 
six-fold degenerate due to manganese spin S = 5/2, 
there are five electrons on the 3d shell and the host 
bands are filled up to the Fermi energy. In the 
excited state, one electron is transferred from the 
host band to the 3d shell (or from the 3d shell to 
the band) resulting in one additional hole (or elec¬ 
tron) in the host bands, respectively. Such virtual 
transfers of electron were used earlier in theories 
of sp — d exchange interaction in semiconductors 
Hill llil- As a consequence of this 3d shell - band 


states hybridization, the ground state of the system 
is no longer six fold degenerate, and the resulting 
low energy excitation spectrum and the anisotropic 
properties of the ion take place, in accord with ex¬ 
perimental observations. This mechanism was de¬ 
scribed in Ref. jl5] and successfully applied to 
EPR results for PbTe:Mn layers grown on BaF 2 and 
KC1 substrates [h|, to the single ion anisotropy in 
Sni-^-yPbyMn^Te mixed crystals 0 and to the 
analysis of magnetic specific heat in Pbi-^Mn^Te 
E3- This method is also used in the present paper 
for calculations of effective spin Hamiltonian. Here, 
we do not repeat detailed description of the method 
and formulas contained in Ref. [15j but sketch only 
the main idea of calculations. 

We consider GeTe crystal at the temperature 
T = 0 with one germanium atom replaced by man¬ 
ganese atom. It is assumed that this replacement 
does not change appreciably neither the band struc¬ 
ture nor the spatial electron distribution. The only 
difference comparing to the pure GeTe is the pres- 
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ence of additional 3d shell on one of the cation 
places with five electrons on it. If we neglect hy¬ 
bridization between 3d orbitals and the orbitals of 
the neighbouring atoms, in the ground state of such 
a system there are five electrons on 3d shell, and 
free electrons fill the band states up to the Fermi 
energy. This ground state is sixfold degenerate due 
to the total spin S = 5/2 of the 3d shell. In the ex¬ 
cited states of the system one electron is transferred 
from 3d shell to an unoccupied band state, or one 
electron from the band is transferred to the 3d shell. 
Such transfers require excitation energy denoted by 
eo, which will be the only fitting parameter for com¬ 
parison of MAE given by the effective Hamiltonian 
with the DFT results. 

In the excited states there are four or six elec¬ 
trons on the 3d shell, and according to the Hund’s 
rule their total angular momentum L = 2 and the 
total spin S = 2. Due to nonzero total angular 
momentum, in the excited states the internal Mn 
spin-orbit coupling ~ L • S should be taken into ac¬ 
count. The inclusion of the internal S-0 coupling 
is crucial for the splitting of the ground state and 
magnetic anisotropy properties of the ion, otherwise 
the ground state is not split, independently whether 
the band spin-orbit coupling is present or not. 

In the Hamiltonian of the system the ground and 
excited states are connected by hybridization be¬ 
tween 3d(Mn) orbitals and orbitals of neighbour¬ 
ing Te ions. The effective spin Hamiltonian for the 
ground state is calculated within the framework of 
the second order perturbation theory for degenerate 
spectra with respect to hybridization. 

The main differences between the present cal¬ 
culations and those presented in Ref. [15] for 
Pbi-^Mn^Te are related to the tight binding model 
of the band structure and to the hybridization pa¬ 
rameters between 3d states of Mn ion and the or¬ 
bitals of neighboring Te ions used in calculations. 
In Ref. [H| we considered Mn ion in PbTe using the 
tight binding model that takes into account near¬ 
est cation - anion neighbors and nearest and next 
nearest cation - cation and anion - anion neigh¬ 
bors. The model was parametrized according to 
Ref. [H]. The integrals describing hybridization Mn 
3d -Te orbitals were taken from Ref. [l9|- Here, 
the tight binding model takes into account inte¬ 
grals up to the third neighbour, and their values 
were taken from DFT calculations for GeTe as the 
matrix elements of Kohn Sham Hamiltonian be¬ 
tween pseudoatomic orbitals of Ge and Te. Using 
ab initio OpenMX package [2 gJ we performed cal¬ 


culations for both fully relativistic and scalar rel¬ 
ativistic pseudopotentials. Thus, we can compare 
effective spin Hamiltonian obtained for situations 
with and without the band spin-orbit interaction. 

From DFT calculations performed for 2x2x2 
supercell containing 64 atoms with one Ge replaced 
by Mn atom we obtain hybridization parameters 
between 3d(Mn) states and s and p orbitals of the 
nearest Te ions. 

In calculations, like in the paper I, the crystal 
structure parameters were assumed ao = 5.98 A, 
a = 88° and r = 0.03. 

Having the effective spin Hamiltonian we ob¬ 
tain Ea^-iP), i- e., the dependence of single ion 
anisotropy energy on the direction of Mn spin de¬ 
fined by the angles 6 and ip. To this aim we intro¬ 
duce a matrix D{9 7 ip) 

D(0,<p) = D z (ip)D y (9), (3) 

where D z and D y are rotation matrices for spin 
S = 5/2 around the z and the y axis, respectively. 
Now, we denote by |a > a state of the Mn spin for 
which the projection on the quantization axis z is 
S z = 5/2, i.e., in this state the Mn spin points in 
the z direction. Then |a' >= D(9,tp)\a > denotes 
the state in which the spin points in the direction 
z' defined by the angles 9 and p. We define 

E A {9,p) =<a'\H\a' > (4) 

as the angular dependence of the single ion 
anisotropy energy on the spin’s direction. For ex¬ 
ample, for the Mn ion in perfect cubic environment 
the effective spin Hamiltonian, Eq. ©, leads to 

E A (9 7 p) = [(cosd) 4 + (sind) 4 ((cos</?) 4 + (siny>) 4 )] 

(5) 

which is identical to that obtained by Yosida and 
Tachiki using a different approach [10]. 

2.2. Determination of transfer energy eg 

To determine the energy eo, i. e. the energy 
necessary to transfer an electron between the 3d 
shell and the top of the valence band, we performed 
scalar relativistic calculations for GeTe. Next, us¬ 
ing the tight binding parameters from the obtained 
Kohn - Sham Hamiltonian we calculated the effec¬ 
tive spin Hamiltonian and the angular dependence 
of MAE, E a {9 7 ip), for for hole concentrations p = 0 
and p = 10 21 cm' 3 for different values of eo- The 
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Figure 1: Angular dependence of MAE for p = 0 (a) and 
p = 10 21 cm -3 (b). Continuous lines correspond to DFT cal¬ 
culations. Dotted and broken lines are single ion anisotropies 
for two values of transfer energy eo- 


Figure 2: Probability density for the angle between easy 
axis and [111] direction (a) and for anisotropy energies (b) 
for different levels of disorder quantified by the spreads of 
the angles 80, 8tp and bond length 8d. 


results were compared to MAE obtained from DFT 
calculations for 64-atom supercell containing one 
Mn ion for vanishing S-0 coupling for Ge and Te 
atoms. Such a comparison makes sense only when 
the band spin-orbit interaction is neglected because, 
as we know from paper I, only in such a case MAE 
is due to Mn ion and does not depend on the direc¬ 
tion of spin polarization of Ge and Te ions. From 
DFT calculations we obtain a 2 C = —0.053 K and 
a 2 C = —0.101 K for insulating and conducting cases, 
respectively. To obtain these values by the effective 
Hamiltonian method we must assume eo=2.92 eV 
for the insulating, and eo=2.58 eV for the conduct¬ 
ing cases, respectively, see the Fig. 03 Because 
in the following we consider the conducting case, 
we assume eo=2.58 eV. The value of eo obtained 
from fitting of the present model to DFT ab initio 
calculations are quite reasonable, and comparable 
with the values used previously in the literature for 
PbTe. In Ref. |ll| eo was considered in the interval 
from 1.6 eV to 3.5 eV, in Ref. [l~3| value 3.5 eV was 
assumed. 

2.3. Dependence of magnetization easy axis direc¬ 
tion on disorder 

To simulate microscopic disorder in Gei-^Mn^Te 
we start from the perfect GeTe lattice character¬ 
ized by do = 5.98 A, a = 88°, r = 0.03 with one 
Ge replaced by Mn. The directions and lengths of 
six bonds connecting Mn with its nearest Te neigh¬ 
bours are specified by (0°, ipf, d°), i = 1,..., 6. We 
assume that the microscopic disorder leads to small 


changes in (0j, <Pi,di), 

0* = 0° + rf <50 

Pi = Pi + rf Sip (6) 

di = + rfSd 

where rf , rf , rf are homogeneously distributed ran¬ 
dom numbers from the interval (-0.5, 0.5). In cal¬ 
culations it is assumed that Mn-Te hybridization 
parameters change with the bond’s length as d ~ 7 / 2 

For a given configuration of Mn surrounding 
we calculate effective spin Hamiltonian, and using 
Eq. Q we find Ea(0,ip). Having this functional 
dependence we obtain direction of the easy axis 
of magnetization and A Ea, the difference between 
maximal and minimal values of Ea(0 , ip). In Fig. [2] 
we present normalized to unity histograms (prob¬ 
ability densities) of directions of the easy axis of 
magnetization with respect to [111] crystallographic 
direction, Fig. 0i, and of AEa, Fig. 0>. These his¬ 
tograms result from 4 x 10 5 random configurations. 

To obtain typical deflections of directions and 
lengths of Mn-Te bonds from those for the ideal 
lattice we performed, using OpenMX package (20l |. 
geometry optimization for the 2 x 2 x 2 64-atom 
supercell. We repeated the optimization for three 
different cases: 1) with one Mn ion, 2) with two Mn 
ions placed at nearest neighbours positions and 3) 
with two Mn ions placed at the two diagonal ver¬ 
tices of the elementary cell and with one Ge vacancy 
on the face of the elementary cell. 

Comparing the directions and lengths of Mn-Te 
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bonds between cases 2) and 3) with those in 1) 
we find that the the deflections are very small. 
The highest deflections are of the order of 1° and 
changes of bond lengths are of the order of 0.01 A. 
We see that even for very small deflections of di¬ 
rections and lengths of bonds from the ideal values, 
the influence on the distribution of the directions 
of the easy axes is quite substantial. However, the 
probability distribution for A Ea for such small 50, 
5<p and 5d remains well centered around the mean 
value, which is of the order of 0.5 K with the vari¬ 
ance less than 0.1 K. 

We have also calculated the averaged over dis¬ 
order the magnetic field dependent magnetizations 
of Mn ion with magnetic field along [111] and per¬ 
pendicular to [111] direction in different tempera¬ 
tures. The differences are noticeable, however in 
very small magnetic fields and extremely low tem¬ 
peratures only. The differences which might be 
compared to those experimentally observed are only 
for unphysically high level of disorder, see Fig. [5] 

We conclude therefore that the single ion 
anisotropy does not play a decisive role in magnetic 
anisotropy energy in Gei_ x Mn x Te. The present 
calculations further confirm the conclusion of the 
paper I that MAE in Gei_ x Mn x Te is mainly due 
to polarization of valence band carriers. 

3. Supercells containing more than one Mn 
ion 

In the present Section we analyze magnetic 
anisotropy for all possible configurations of two Mn 
ions in the 2 x 2 x 2 supercell, for the case of 
three Mn ions we show only results for a few ran¬ 
domly chosen configurations. Technical details of 
calculations were presented in the paper I. We used 
OpenMX package [20]. In all calculations of MAE 
the number of integration points in the Brillouin 
zone were 512 and 2744 for insulating and conduct¬ 
ing cases, respectively. As was shown in the paper 
I such numbers give convergent results. In the con¬ 
ducting case, the calculations were performed for 
hole concentration p = 10 21 cm~ 3 , which is typical 
for samples studied experimentally. For geometry 
optimization, the number of integration points in 
the Brillouin zone reduced to 64 give convergent 
results. The force criterion for geometry optimiza¬ 
tion was equal to 5 x 10 ~ 3 Hartree/Bolir. 

In order to describe possible configurations, first, 
we introduce notation describing positions of Mn 
ions in the supercell. 


Let us consider three vectors Ci, C 2 and C 3 , parallel 
to the [ 100 ], [ 010 ] and [ 001 ] crystallographic direc¬ 
tions, respectively. The lengths of the vectors are 
the same, ao/2. The symbol {ijk} describes the po¬ 
sition of an atom in the supercell r = ici +jc 2 + kca 
where i, j, k are integer numbers. 

The first Mn atom is always placed at {000}. The 
2x2x2 supercell contain 32 different cation po¬ 
sitions, thus the second Mn atom may be placed 
at 31 remaining positions. However, some of the 
resulting 31 configurations are equivalent because 
of the periodic boundary conditions; for example 
configurations {110} and {330} are equivalent. Ac¬ 
tually, there are 21 nonequivalent configurations. 
We stress that due to nonzero Mn spin this num¬ 
ber is larger then the number resulting from purely 
geometrical symmetry considerations for a rhombo- 
hedral lattice. Consider, for example, three nearest 
neighbor pair configurations for which the second 
Mn atom is placed at {110}, {011} and {101}, re¬ 
spectively. For vanishing spin, those three configu¬ 
rations are equivalent from the purely geometrical 
point of view, being related by the rotation along 
[111] axis by the angle ±120°. However, for finite 
S, they are nonequivalent if we consider magnetic 
properties of the crystal, what is explicitly shown 
in the following in Fig. [3] 

In the calculations, the Mn spin vectors are 
placed in the (-110) plane. Due to reflection sym¬ 
metry with respect to this plane, the configurations 
{ 011 } and { 101 } give identical results for the angu¬ 
lar dependence of E(0), thus it is enough to perform 
calculations for one of them, but for arbitrary direc¬ 
tion of manganese spins it is necessary to calculate 
MAE for both. 


In Figs. [3]and|3]we present the results for the cho¬ 
sen example configurations for insulating and con¬ 
ducting cases, respectively. Also we show arithmeti¬ 
cal averages over all nonequivalent configurations. 
Let us notice, however, that the results presented 
in the paper I indicate that even in the insulating 
case the pair of Mn ions in the 2x2x2 supercell 
cannot be treated as isolated. Thus, strictly speak¬ 
ing, these arithmetical averages cannot be treated 
as the averages over the quenched disorder. Never¬ 
theless, they give some information about behavior 
of MAE in macroscopic crystal where the magnetic 
ions are placed randomly. 

In Figs. [3] and 0] we compare also MAE for nonre- 
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Figure 3: Angular dependence of the energy of magnetic 
anisotropy for 2x2x2 supercell containing two Mn ions 
for insulating case. Continuous lines - nonrelaxed lattice, 
broken lines - relaxed lattice. 



Figure 5: Angular dependence of MAE for several config¬ 
urations of three Mn ions in 2 x 2 x 2 supercell for hole 
concentration p = 10 21 cm -3 . 



Figure 4: Angular dependence of the energy of magnetic 
anisotropy for 2x2x2 supercell containing two Mn ions 
for hole concentration p = 10 21 cm -3 . Continuous lines - 
nonrelaxed lattice, broken lines - relaxed lattice character¬ 
ized by lattice parameters ao = 5.98 A, a = 88°, r = 0.03. 
The dotted lines are for nonrelaxed lattice for other lattice 
parameters ao = 5.98 A, a = 88.3°, r = 0.025. 


Table 1: Values of coefficients a^ c and a 2 s describing MAE 
for 2x2x2 supercell containing two Mn ions for nonrelaxed 
and relaxed lattices and hole concentration p = 0. 




nonrelaxed 

relaxed 

configuration 

/ 

fl2c 

«2s 

«2c 

02 s 

{ 110 } 

i 

0.49 

0.69 

0.62 

0.92 

{ 011 } 

2 

0.48 

-0.34 

0.20 

-0.47 

{ 002 } 

1 

- 0.12 

-1.13 

-0.84 

-1.56 

{ 112 } 

1 

0.24 

- 0.01 

- 0.11 

- 0.10 

{013} 

2 

-0.33 

-0.42 

-0.72 

-0.36 

{ 020 } 

2 

0.31 

0.57 

0.42 

0.80 

{130} 

2 

0.04 

0.96 

0.05 

0.71 

{ 121 } 

2 

0.22 

0.01 

-0.05 

0.06 

{ 022 } 

2 

0.32 

-0.16 

- 0.10 

-0.27 

{132} 

2 

0.36 

- 0.22 

-0.05 

- 0.21 

{123} 

2 

0.38 

0.11 

0.15 

0.11 

{ 220 } 

1 

0.26 

0.29 

0.15 

0.50 

{ 222 } 

1 

0.48 

- 0.01 

0.20 

0.00 

average 


0.23 

0.04 

- 0.01 

0.02 
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Table 2: Values of coefficients a 2c and a 2s describing MAE 
for 2x2x2 supercell containing two Mn ions for nonrelaxed 
and relaxed lattices and hole concentration p = 10 21 crn . 




nonrelaxed 

relaxed 

configuration 

/ 

«2c 

d2s 

U2c 

U2s 

{ 110 } 

i 

-3.91 

4.6 

1.88 

1.03 

{ 011 } 

2 

-6.93 

-2.51 

-3.02 

-0.61 

{ 002 } 

1 

-3.81 

0.68 

-1.55 

2.14 

{ 112 } 

1 

-7.31 

-1.92 

-6.39 

- 1.12 

{013} 

2 

-3.02 

-1.84 

-4.22 

-2.96 

{ 020 } 

2 

-3.13 

-0.46 

1.27 

- 1.02 

{130} 

2 

-0.29 

3.74 

-0.42 

5.92 

{ 121 } 

2 

-7.16 

1.09 

-7.16 

0.98 

{ 022 } 

2 

-7.47 

-3.22 

-6.56 

-3.96 

{132} 

2 

-6.51 

-2.35 

-4.08 

-2.16 

{123} 

2 

-5.67 

1.06 

-2.42 

1.53 

{ 220 } 

1 

-4.87 

6.39 

-2.84 

8.11 

{ 222 } 

1 

-3.36 

0.11 

-1.48 

0.48 

average 


-4.93 

0.04 

-3.03 

0.29 



Figure 6: Angular dependence of contributions of kinetic 
(a), electrostatic (b) and exchange-correlation (c) energies 
to magnetic anisotropy energy (d) for 2x2x2 supercell 
containing two Mn ions in {110} (continuous lines) and 
{130} (broken lines) configurations. The hole concentration 
is p = 10 21 cm -3 . 


laxed (continuous lines) and relaxed (broken lines) 
supercells. As in the case of single Mn ion in the su¬ 
percell, the relaxation of the lattice leads to signifi¬ 
cant changes in MAE. For comparison, by the dot¬ 
ted lines we show MAE for nonrelaxed lattice with 
slightly different lattice parameters: ao = 5.98 A, 
a = 88.3°, t = 0.025. Like in paper I, for single Mn 
ion in the supercell, such changes result in changes 
of MAE. 

The results for all configurations are summarized 
in Tables [T| and [5] The first column of each table 
gives configurations, in the second column there are 
corresponding weight factors /. The factor / = 1 
is for a configuration which remains the same when 
reflected in the (- 110 ) plane, and / = 2 if there exist 
two nonequivalent configurations related by reflec¬ 
tion in the (-110) plane. The coefficients a 2 C and 
ct 2s describe curves used to fitting the calculated 
points: 

Ea(0) = a 2c cos (2 6) + a 2s sin(20) (7) 

It turns out that inclusion of terms proportional to 
cos(40) and sin(40) to the fitting formula does not 
improve quality of the fit, because the coefficients 
describing these terms are more than order of mag¬ 
nitude smaller than a 2c and a 2s . 

From the presented results we draw two main 
conclusions. 

First, for a given pair configuration the direction 
of the easy axis of magnetization, in general, is 
very different from [ 111 ] crystallographic direction 


or the direction perpendicular to [111]. However, 
after the assumed averaging procedure the average 
of the a 2s coefficient is very small, it means that 
the average MAE possesses rhombohedral symme¬ 
try of the macroscopic crystal. Similar behavior 
was previously observed in calculations of single ion 
magnetic anisotropy in Sni-^Mn^Te Q in which 
disorder effects were taken into account. 

Second, in most of cases the values of yj a| c + a| s 
for pair configurations are significantly larger than 
the average of coefficient a 2c . This fact is also con¬ 
nected with varying directions of the easy axes of 
magnetization for different pair configurations. 

In Fig. [5] we show results for several randomly 
chosen configurations of three Mn ions in 2 x 2 x 2 
supercell. We are not able to perform calculations 
for all possible configurations because in this case 
the number of configurations is too large, however 
one may notice that the qualitative results are sim¬ 
ilar to the case of two Mn ions in the supercell. 

Angular dependencies for the kinetic, electro¬ 
static and exchange-correlation components of 
MAE are shown in Fig. [ 6 ] Their nonconstant an¬ 
gular dependencies prove that the spatial charge 
distribution varies with the Mn spin’s direction. 

Two features may be noticed. First, the changes 
of kinetic, electrostatic and exchange-correlation 
energies with the direction of manganese spins are 
much larger than the changes of their sum, the mag¬ 
netic anisotropy energy. However this is usual in 
DFT calculations - the dependence of components 
















of the total energy on external parameters is much 
stronger than the dependence of the total energy 
itself. The other feature is more interesting. Let 
us notice that the maxima and minima of compo¬ 
nents are placed in different positions than the max¬ 
ima and minima of MAE. This was not the case 
for single manganese ion in a supercell considered 
in the paper I where the maxima or minima for all 
curves were at the same positions. This shows, even 
more expressively than in the paper I, the close re¬ 
lationship between spin direction of manganese ions 
and the spatial distribution of electron charge. The 
changes of spatial charge density are caused by the 
spin-orbit interaction. 

4. Conclusions 

In this work, we theoretically analyzed the influ¬ 
ence of chemical disorder on magnetic anisotropy 
in ferromagnetic Gei-^Mn^Te semiconductor lay¬ 
ers known to undergo ferroelectric structural dis¬ 
tortion of rhombohedral crystal lattice and exhibit 
perpendicular magnetic anisotropy. Using DFT 
method we determined the local crystal structure 
in Gei-^Mn^Te substitutional solid solution and 
calculated the single ion contribution to magnetic 
anisotropy energy. Our calculation revealed that 
the single ion contribution is order of magnitude 
smaller than the contribution originating from spin 
polarization induced by Mn ions into neighboring 
Te and Ge ions. This effect results from the hy¬ 
bridization of magnetic 3 d orbitals of Mn ions and 
valence band states of GeTe subject to spin-orbit 
interactions Q. 

We also discussed magnetic anisotropy effects due 
to pairs of Mn ions differently allocated in the su¬ 
percell and showed that the spatial averaging over 
the chemical disorder strongly reduces the strength 
of this magnetic anisotropy mechanism and restores 
the global rhombohedral symmetry of magnetic sys¬ 
tem. 

Finally, the calculations show also that although 
in principle ab initio numerical methods may be ap¬ 
plied to the calculations of MAE, in order to take 
into account disorder and make a direct quantita¬ 
tive comparison with experiment, much larger sys¬ 
tems must be considered. 
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